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Abstract 

With multi-stream coannular nozzles, sometimes tones occur that may cause the nozzle to fail noise 
regulation standards. A two-stream nozzle was studied experimentally and numerically in an at- 
tempt to identify the sources of such tones and explore remedies. For the given nozzle configura- 
tion, sharp tones occurred in a range of low jet Mach numbers. The tones apparently occurred due to 
a coupling between vortex shedding from the struts, which held the nozzles and the center-body to- 
gether, with various duct acoustic modes. A leading edge treatment of the struts is shown to elimi- 
nate the tones via disruption of the vortex shedding. 


Introduction 

There have been numerous previous studies addressing tones occurring with jets and shear layers 
and their interactions with structures [e.g., 1-5]. The present investigation addresses tones some- 
times observed with multi-stream nozzle flows. In fact, the study was prompted by persistent tones 
encountered in large-scale multi-stream nozzle tests at NASA Glenn Research Center (GRC) [6]. 
Such tones may make the nozzle fail noise regulation standards and this led to a model-scale exper- 
imental investigation. A 0.46-scale model of a two-stream coannular nozzle, that exhibited such 
tones in the large-scale tests, was fabricated for the investigation. The effort was to understand the 
origin of such tones and explore possible remedies. This paper reports results from the model-scale 
experiment. 


The flow geometries of complex multi-stream nozzles involve multiple coannular ducts as well as a 
center-body. Usually, three or four struts, spaced equally in the circumferential direction, hold these 
coannular elements together. As the flow velocity is increased from zero, distinct tones may occur. 
With increasing speed a tone may disappear and new tones with different frequencies may appear. 
As found in this investigation and as explained with the results in the following, one mechanism for 
the origin of such tones is vortex shedding from the struts. The shedding apparently couples with 
acoustic resonance modes of the nozzle configuration. The result is an acoustic tone that 1s spectral- 
ly more concentrated than that from vortex shedding alone. A remedy is found. A leading edge 
treatment of the struts disrupts organized vortex shedding via a disruption of the two-dimensionality 
of the flow. This may yield a complete elimination of the tones. Key results from the investigation 
are summarized in this paper. 
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Experimental Facility 

The experiments are conducted in an open jet facility at GRC. A picture of the nozzle mounted on 
the end of the facility is shown in Fig. 1(a). Unless specifically stated all dimensions are given in 
inches. A 30" diameter plenum chamber with flow conditioning units 1s the main component of the 
facility. Compressed air supplied from one end of the plenum discharges through the nozzle, 
mounted on the other end, into the ambient of the test chamber. The two-stream nozzle assembly is 
shown by the drawing in Fig. 1(b), with annotations identifying the various components. It consists 
of a cylindrical inlet section that is attached to the plenum chamber. This section 1s of 4.66" diame- 
ter and 3.5" in length. The diameter and slope at the junction with the plenum are matched for 
smooth flow lines. An interested reader may find further description in [7]. 


With reference to Fig. 1(b), the rest of the two-stream nozzle assembly is bolted to the end of the 
cylindrical inlet section. The outer nozzle has an exit diameter of 3.46". Compressed air from the 
single supply splits and passes through the inner and outer passages. The inner flow passes through 
a removable nozzle having exit diameter of 3.03". Three otherwise identical inner nozzles are used 
whose lip thickness (typ) 1s varied as 0.015", 0.030" and 0.050". Items 1 and 2 in Fig. 1(c) represent 
the inner nozzle with the thinnest and the thickest lip thicknesses, respectively. The ‘removable in- 
let’ in Fig. 1b is a straight cylindrical section with rounded leading edge. Most of the data pertain to 
an inlet with length L= 2". Two other inlets of length L= 0.75" and 4.75" (items 3 and 4 in Fig. Ic) 
are also used in order to vary the (longitudinal) duct acoustic modes and study their impact on the 
audible tones. The ratio of the flow velocities at the exit of the inner and outer channel is measured 
to be approximately unity for all hardware at all operating pressures. Items 5 and 6 in Fig. 1(c) are 
two annular discs that can be placed just upstream of the struts to block either the inner flow or the 
outer flow. 


The inner nozzle and the center-body are held together with respect to the outer shell of the assem- 
bly by four struts. The struts are placed in a ‘crossed’ shape. Thus, both the outer and the inner 
flows encounter four struts spaced equally in the circumferential direction (Fig. 1b). A sketch of a 
strut spanning the inner channel is shown in Fig. 1(d). The span (s) is 0.8" and the chord (c) 1s 
0.625" while the maximum thickness (f) is 0.125". Also shown in this sketch is a leading edge 
‘treatment’ in the form of a cap attached to the strut at mid-span. The cap geometry is described fur- 
ther with the results in the following. Results from numerical analysis of acoustic waves within the 
nozzle, and associated procedure, are discussed in a separate subsection following the experimental 
results. 


Experimental Results 

Sound pressure level (SPL) spectra measured at a location perpendicular to the nozzle axis, and 
about 50" away from the exit, are shown in Fig. 2. The spectral traces (O-100kHz range with 1000- 
line analysis) are staggered by 10 dB and the jet Mach number (M)) is indicated for each trace. A 
range of Mj from low subsonic to just over unity is covered. These data are for an inner nozzle (Fig. 
1b) with lip thickness of 0.030". The spectra display a broadband peak; e.g., at about 25 kHz at M; = 
0.256 and 90 kHz at M; =1.116. Its frequency continuously increases with increasing Mj. These 
peaks are due to vortex shedding from the lip of the inner nozzle, as evident from corresponding 
non-dimensional value of the frequency. The frequency (fs) of the peak is plotted in Fig. 3 as Strou- 
hal number based on the lip thickness and jet velocity. Also shown in this figure are corresponding 
data for two other inner nozzles with lip thicknesses of 0.015" and 0.050". The Strouhal number 
turns out to be about 0.20, a value that corresponds to the non-dimensional frequency of von 


Karman vortex shedding. The rest of the data in this paper pertain to the inner nozzle with 0.030" 
lip thickness. 


In addition to the broadband peak due to vortex shedding from the lip, sharp spikes are noted in the 
spectra at lower Mach numbers (Fig. 2). These are accompanied by audible tones. The trends of 
those tones are shown more clearly in Fig. 4 for the same nozzle configuration with spectral traces 
measured with a narrower frequency range (O-20kHz with 1000-line analysis). As the jet Mach 
number is increased the frequencies of the tone vary in stages. That is, the frequency remains con- 
stant over a range of M; then jumps to a new value and the pattern continues. This is further demon- 
strated by the tone frequency data in Fig. 5. There are four major stages as marked in the figure. 
With the calculated velocity within the inner channel near the struts (~0.388/M;) and the strut thick- 
ness, t=0.125" (Fig. 1d), the dimensionless frequency (Strouhal number) in each stage turns out to 
be about 0.3. This is not too far from the value for Karman vortex shedding. (The Strouhal number 
is often found over a range around the expected value of about 0.2. This is due to the fact that the 
effective length scale is different from what is used as well as boundary layer effects. In the present 
case, a thickness smaller than ¢ corresponding to a location closer to the trailing edge is likely to be 
the length scale, accounting for the higher value of 0.3 that was obtained by using f). The reader is 
reminded that, unlike the broadband peaks in Fig. 2 due to shedding from the inner nozzle lip, here 
we are considering shedding from the trailing edge of the struts. The sharp peaks at lower Mj, corre- 
sponding to the audible tones, are the focus of the present study. 


In an effort to understand the nature of the acoustic modes, spectral data were acquired with varying 
length of the inlet section (Fig. 1b). Corresponding frequency data for the tones are plotted in Fig. 6. 
Whereas the data in Fig. 5 corresponded to an inlet length L=2", data for L=0.75" (short) and 4.75" 
(long) inlets are shown in Fig. 6. Also included are data with the outer flow blocked using the annu- 
lar disc 6 of Fig. I(c). It is clear that for any of these parametric variations the tone frequencies and 
the staging behavior basically remain unchanged. When the inner flow was blocked by disc 5, Fig. 
l(c), the tones disappeared. Thus, the source of the tone must be in the inner flow since the tones 
occur regardless of the presence of the outer flow. Furthermore, somewhat surprisingly, it becomes 
clear that the length of the inner duct has no impact on the frequencies or the stages of the tones. 
The duct modes are addressed further with the numerical results in the following. All subsequent 
data in the paper pertain to the inlet section with length L=2". 


From the discussion of the data in Figs. 4-6 it is apparent that vortex shedding from the struts is the 
likely origin of the observed tones. However, the fact that the tone frequency remains constant with- 
in each stage (Figs. 5, 6) indicates that there is a coupling with acoustic resonant frequencies char- 
acteristic of the nozzle configuration. When the shedding frequency is close, a particular acoustic 
resonance mode 1s excited. The excitation of that mode occurs over a range of the jet Mach number. 
When the velocity is varied sufficiently the shedding frequency becomes incompatible and the cou- 
pling ceases. However, the frequency can now become compatible with another acoustic mode lead- 
ing to a tone of different frequency. This explains the staging behavior. 


It should be mentioned here that various other parametric variations were examined in an effort to 
understand the origin of the sharp tones. The parameters varied included the lip thickness of the in- 
ner nozzle, lip-to-lip distance between the inner and the outer nozzle, a flared and a constricted inlet 
instead of the straight removable inlet in Fig. 1b, as well as the overall length of the inner duct by 
varying the length of the removable inlet as already discussed. (When designing the experiment, the 


origin of the tones was suspected to be with the flow near the nozzle exit and its interaction with the 
duct modes. The struts were not even considered as culprits and hence no provision was made for 
varying the strut geometry.) For these explorations four conditions were chosen from the four stages 
of Fig. 5. Corresponding spectral traces at M; = 0.168, 0.260, 0.345 and 0.427 are shown in Fig. 7, 
which display spikes at frequencies 4.5, 7.75, 9.5 and 12.38 kHz, respectively. For clarity the spikes 
are marked by the circular symbols. Variation in the lip-to-lip distance was achieved by unscrewing 
the inner nozzle (that, of course, left a cavity on the inner wall of the inner duct). None of these var- 
lations affected the two spectral peaks at the lowest frequencies (4.5 and 7.75 kHz, Fig. 7). With lip- 
to-lip distance variation the spectral peaks at higher M; were altered with appearance of peaks at 
other frequencies, probably due to additional tones originating from the inner wall cavity. In gen- 
eral, the tones remained unchanged for variation of the aforementioned parameters and eventually it 
became apparent that the tones were tied to the struts in the inner duct. 


A key result of the study was obtained when the leading edge cap was applied on the inner struts. 
The tones were completely eliminated as demonstrated in Fig. 8. The caps were cut from a plastic 
‘edge protector’ for car doors. With reference to Fig. 1(d), the dimensions were approximately, /= 
0.31", h=0.26" and w=0.64". These were placed centrally on each of the four struts in the inner 
flow. Experiments showed that placing the caps on the struts in the outer annular flow had no effect 
on the tones. The tones were suppressed only when the caps were placed on the inner struts. Fur- 
thermore, it was found that caps with much smaller span w= 0.3" or even 0.1" were just as effective 
in eliminating the tones. Apparently, the caps disrupted the two-dimensionality of the flow over the 
struts and thereby the organized vortex shedding was stopped which in turn led to the elimination of 
the tones. 


A further interesting result was obtained when caps with full span, equal to the span of the struts 
(w= 0.8"), were tried. With the full-span caps on all four inner struts the tones reappeared albeit at 
higher Mach numbers. This is shown in Fig. 9. A tone at 4.13 kHz appeared at the highest Mj. At 
even higher values of M; two stages of the tone could be detected, as can be seen from the spectra 
traces of Fig. 10. With the no-cap condition tones were not heard at M; values higher than about 
0.43. With the full-caps the occurrence of the tones shifted to a higher Mj; range and could be heard 
up to about Mj= 0.85. With the height h=0.26" of the cap (Fig. Id) and the local velocity in the duct 
(~0.388M;), the frequencies of the two stages are found to correspond to a Strouhal number of about 
0.22. Therefore, it is apparent that with the full-span caps two-dimensionality of the flow is restored 
that leads to organized vortex shedding which in turn excites the tones. However, since the shedding 
frequencies are lower due to the larger thickness, for a given acoustic mode the frequency matching 
occurs at a higher Mach number. This is why the occurrence of the tones has shifted to a range of 
higher Mj. 


For the full-caps the frequency data are plotted in Fig. 11. Also shown in this figure are the ampli- 
tudes of the spectral peaks (in arbitrary scale) to provide an idea about the relative variation of the 
tone intensity within each stage. Clearly, the amplitude becomes the largest in the middle of each 
stage. Also, these amplitudes at higher M; with the full-caps turned out to be much larger than those 
encountered with the no-cap case. 


The higher Mach number conditions with the full caps allowed good visualization of the flow using 
the schlieren technique. Typical visualization pictures for the two stages of Fig. 11 are shown in 
Fig. 12. The large amplitude tones clearly excite the outer shear layer as evident from the organized 


vortices marked by the downward pointing arrows. Vortex shedding from the lip of the inner nozzle 
(leading to the broadband peaks, Fig. 2) can be barely discerned as marked by the upward pointing 
arrows. Unfortunately, neither these pictures nor the experimental results presented so far provide a 
clear idea about the nature of the acoustic modes within the nozzle that are excited by vortex shed- 
ding from the struts. Some insight is gained by the numerical simulation results described next. 


Numerical simulation 

The intent of this analysis is to identify duct modes that couple with Strouhal shedding from the 
struts. A commercial software package, COMSOL Multiphysics® [8], 1s used to perform the analy- 
sis for the 2-stream nozzle configuration. A Computer Aided Design (CAD) geometry (Fig. 1b) 1s 
used to create a 3D COMSOL model of the nozzle. The computed domain (Fig. 13) starts upstream 
of the internal flow split between the core- and bypass-flow paths, and extends downstream to just 
past the tip of the center plug. The simulation solves the acoustic pressure physics using a finite el- 
ement formulation on an unstructured mesh. There is no flow in the domain. The upstream bounda- 
ry has a plane wave radiation condition. The downstream radial and axial surfaces use a spherical 
wave radiation boundary condition with center on the jet axis and axially at the inner nozzle lip. 


Two small panels on either side of one of the struts in the inner flow are used to provide periodic 
perturbation. The panels are placed near the trailing edge. The “driven strut’ at the top, with the 
panels marked blue, is identified in Fig. 13. The two sides are driven 180 degrees out of phase in 
order to simulate the asymmetric vortex shedding. The perturbation pressure amplitude is set arbi- 
trarily at + 1 Pascal. 


With the perturbation at a given frequency the simulation is performed and the acoustic field in the 
domain is then interrogated. For each frequency, two metrics for resonance strength are collected: 
the maximum of the pressure fluctuation amplitude in the domain and the maximum of the root 
mean square acoustic velocity in the domain. A sweep in frequency is then conducted in steps, e.g., 
25Hz, 100Hz, etc., with larger steps being used at higher perturbation frequencies. This way a spec- 
trum of the response function is constructed that reveals the resonance conditions. Then, for specific 
resonant frequencies the mode shape, 1.e., the nature of the pressure distribution in the domain is 
studied. As a first step, the mode shape over a cross-sectional plane just downstream of the strut 
(marked red and identified in Fig. 13) 1s examined. 


The response functions (maximum pressure amplitude within the domain and corresponding maxi- 
mum r.m.s. velocity) versus the perturbation frequency are shown in Fig. 14. It immediately be- 
comes clear that the simulation approximately captures the tone frequencies as noted in the experi- 
ment. The three peaks seen in Fig. 14(b) roughly correspond to three of the peaks seen in Fig. 7 at 
4.46, 7.76 and 12.37 kHz. The peak at 9.76 kHz (Fig. 7) is missed, however, there is a hint of addi- 
tional energy around that frequency in Fig. 14(b). Furthermore, the identified resonant frequencies 
are encouraging in that they exhibit amplification at the non-driven struts, indicating a synchronized 
Strouhal shedding. As discussed further in the following, the frequency agreement with the experi- 
ment is not exact, however, quite satisfactory considering the facts that the simulation 1s for no-flow 
and there are experimental uncertainties (approximately 100 Hz) due to variations in ambient condi- 
tions from day to day. 


The pressure distributions on the monitored plane (Fig. 13) are shown in Fig. 15 for the three domi- 
nant resonant frequencies of Fig. 14(b). These frequencies are 4.525, 8.05 and 12.52 kHz in Figs. 
15(a), (b) and (c), respectively. Each distribution represents the amplitude in phase with the driven 
surfaces. Clearly, ‘standing wave’ patterns have set up among the four struts. A region of either 
negative or positive pressure exists between a given pair of struts in Fig. 15(a) and the pattern is an- 
ti-symmetric across a diameter containing the opposite struts. It appears that the azimuthal spacing 
of the struts plays a role in dictating the resonance condition. The pattern in Fig. 15(a) corresponds 
to the tone at the lowest frequency stage (Fig. 5) and thus may be considered as the ‘fundamental’ in 
the phenomenon. 


The data for resonance at 8.05 kHz (Fig. 15b) exhibit pairs of positive and negative pressure regions 
within an intra-strut space. Thus, this may be considered as the first harmonic in the phenomenon. 
The distribution in Fig. 15(c) is also similar but there is an additional radial distribution. This repre- 
sents a complex higher harmonic. Resonances at many such harmonics are possible. In fact, an 
eigenmode analysis using the COMSOL solver identified great many eigenfrequencies. For exam- 
ple, 98 modes were found for the given geometry in the 3500-5500 Hz band. Only a selected few of 
these would be compatible with, and thus excitable by, the asymmetric perturbation from the TE of 
one of the inner struts. Note that due to prohibitive computational time at high frequencies only one- 
half of the cross-section 1s computed for the case in Fig. 15 (c). Furthermore, the amplified scale in 
the latter case (1.5 Pa, compared to 3.5 Pa in Fig. 15b) made the pressure distribution within the 
outer annular passage also visible. 


A corresponding pressure distribution for the fundamental (4.525 kHz) over the entire domain is 
shown by six 1so-surfaces in Fig. 16(a), for pressures of +1, +2 and +3 Pa. Another perspective of 
the same domain is shown in Fig. 16(b), skipping the +2 Pa 1so-surfaces to avoid clutter. These dis- 
tributions indicate a complex mode shape. It 1s not a simple longitudinal, helical or radial duct 
mode. From the +3 Pa iso-surfaces, marked in Fig. 16(b), it is apparent that regions of high pressure 
have formed, against the outer wall of the inner channel, in between pairs of struts. (The walls of the 
inner and outer channels are made transparent in the figures in order to make the iso-surfaces visi- 
ble.) The +1 Pa iso-surfaces (with lighter red and blue colors, Fig. 16b) touch the ‘driven’ strut 
since the driving pressure is of that amplitude. Furthermore, the latter iso-surface touches all the 
struts demonstrating the synchronized motion in the domain. There must be acoustic feedback 
among the struts in order to achieve the synchronization of the vortex shedding. Note once again 
that there is high pressure buildup in the intra-strut spaces. These are locations of pressure anti- 
nodes. Therefore, the struts themselves should be locations of pressure nodes and hence anti-nodes 
for induced velocity. This is indeed found to be the case as demonstrated in Fig. 16(c). High in- 
duced velocity near the TE of the struts would be consistent with synchronized vortex shedding in 
the case with flow. 


The high pressure buildup against the inner path outer wall in the intra-strut spaces also indicate re- 
inforcement by reflected waves from the wall, and hence a necessary role of the duct surrounding 
the struts to establish the standing waves. In other words, such resonance would not be expected if 
the walls were removed and the strut configuration were left in free space. Except for imposition of 
a cut-on frequency for the waves to emanate through the nozzle exit, any additional role of the ducts 
(such as, possible feedback between the struts and the nozzle exit) in establishing the resonances 
remains unclear at this time. In any case, the standing wave pattern around the struts under the arti- 
ficial perturbation appears to be the source from which waves propagate in both upstream and 


downstream directions. In practice, with the flow on, the organized vortex shedding substitutes for 
the artificial perturbation to generate the tones. 


Finally, it should be noted that the resonance frequencies found from the simulation were somewhat 
higher than those seen in the experiment. It is possible this is due to the fact that the simulation in- 
volved no-flow. With flow, the feedback from one strut to the next would require longer wave- 
length, directed upstream at an angle, whose component on the plane of the struts would excite and 
synchronize with shedding from the adjacent struts. Thus, a simple geometric correction to account 
for the longer wavelength due to local flow velocity would yield the correct frequencies. Table 1 
shows such a correction applied to the predicted frequencies from the no-flow simulation and the 
resulting agreement with the experimental data. 


Table 1 Comparison of predicted frequency with tones observed in the experiment. 


M,, Jet Moiocal, Ner «Correction, Freq, f(Hz) Freq, f (Hz) Freq, f | Difference 


Mach num- — channel M Sqrt(1- Mio- Predicted Corrected (Hz) from 
ber at strut loca- sat) Expt Expt. 
tion 
0.168 0.065 0.998 4575 4565 4460 2.4% 
0.260 0.101 0.995 8095 8054 7760 3.8% 
0.427 0.166 0.987 12690 12522 12395 1.2% 
Conclusions 


A two-stream nozzle was studied experimentally and numerically in an attempt to identify the 
sources of resonant tones. With the given nozzle configuration, sharp tones occurred in a range of 
low jet Mach numbers. These appeared due to a coupling between vortex shedding from the struts, 
which held the nozzles and the center-body together, and acoustic resonance modes characteristic of 
the nozzle configuration. With the present nozzle geometry the source of the tone is determined to 
be the four struts within the inner channel. A standing wave pattern occurred around those struts 
whose shape varied and became more complex with increasing frequency. The ‘fundamental’ in- 
volved alternate positive and negative pressure cells in the intra-strut spaces, the pattern being anti- 
symmetric about a diametral plane containing the opposite struts. At higher frequencies, more com- 
plex mode shapes occurred involving multiple pressure cells in a given intra-strut space as well as 
multiple cells occurring in the radial direction. A leading edge treatment of the struts in the inner 
channel is shown to eliminate the tones. This occurred due to a disruption of the organized vortex 
shedding. Apart from the nozzle flows considered in this study, conceivably a similar leading edge 
treatment could be applicable in many other situations, e.g., in wind tunnel tests, where tones in- 
duced by similar vortex shedding might be a problem. The remedy is simple and can be included in 
the original design of the hardware or retrofitted with existing hardware. 
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Fig. 3 Strouhal number corresponding to the 

broadband spectral peaks (Fig. 2) versus Mj, 

for three different lip thicknesses of the inner 
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Fig. 1 Experimental facility. (a) Jet rig, (b) 
drawing of nozzle assembly, (c) accessories, 
(d) inner flow strut with leading edge ‘cap’. 
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Fig. 4 Sound pressure level spectra at differ- Fig. 6 Tone frequency versus M; with inlets of 
ent jet Mach number for the 0.030" lip inner two different lengths and with a case where 
nozzle, in a low M; -range with expanded ab- the outer flow is blocked. 

Scissa. 


Mj=0.427 


et 
le eel * 
Pa whe 


* 
a ga 
ees oe a 


fy (KHz) 
SPL, dB (staggered) 


Fig. 7 Sound pressure level spectra at four jet 


Fig. 5 Tone frequency versus Mj for the nozzle Mach numbers for detailed parametric study. 
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Fig. 10 Sound pressure level spectra with full- 
span caps on the struts in the inner flow, at 
higher values of Mj. 


Fig. 8 Sound pressure level spectra at same 
four M; of Fig. 7 but with semi-span caps on 
the struts in the inner flow. 
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Fig. 9 Sound pressure level spectra at same Fig. 11 Tone frequency versus M; with full- 
four M; of Fig. 7 with full-span caps on the span caps on inner struts. 
struts in the inner flow. 
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Fig. 12 Schlieren pictures of the flow; full- 
span caps on four inner struts. (a) Mj =0.45, 
f=4.13 kHz, (b) Mj =0.67, f=6 kHz. 
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Fig. 13 Computational domain. Monitored 
plane at the trailing edge of struts 1s shown in 
red. Also identified is the driven strut. The 
coordinate labels in this figure are in meters. 
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Fig. 14 Maximum of response function within 
computed domain versus frequency of pertur- 
bation. (a) Maximum pressure, (b) maximum 
r.m.s. velocity. 


Fig. 15 Cross-sectional distribution of pres- 
sure at the monitored plane (Fig. 13). Fre- 
quency and amplitude range are: (a) 4.525 
kHz, +2.03 Pa, (b) 8.05 kHz, +3.5 Pa, and (c) 
12.52 kHz, +1.5 Pa. 
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Fig. 16 Iso-surfaces of pressure and velocity for 
f=4.525 kHz over entire domain. (a) Side view 
with iso-surfaces of pressure at +1, +2 and +3 
Pa. (b) Oblique view with iso-surfaces of +1 and 
+3 Pa. (c) Oblique view with an iso-surface of 
r.m.s. velocity at 0.006 m/s. 


